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Abstract. We extend our result of Q] and show that one can as- 
sociate with the stochastically perturbed non-viscid Burgers equa- 
tion a system of viscous balance laws. The Cauchy data for the 
Burgers equation generates the data for this system. Till the mo- 
ment of the shock formation in the solution to the Burgers equa- 
tion the above system of viscous balance laws can be reduced to 
the pressureless gas dynamics system (in a limit as the parameters 
of perturbation tend to zero). If the solution to the Burgers equa- 
tion contains shocks, the limit system is equivalent to the system 
with a specific pressure, in some sense analogous to the pressure of 
barotropic monoatomic gas. 

1. Stochastic perturbation of the Burgers 
equation 

Let us consider the Cauchy problem for the non-viscous Burgers 
equation: 

(1.1) d t u+(u, V)u = 0, t > 0, u(x, 0) = u (x) e C\R n )nC b (^L n ). 

Here u(x,t) = (ui, ...,u n )(x,t) is a vector-function M n+1 — > M. n . 

We associate with (II. ip the following system of stochastic differential 
equations: 
(1.2) 

dX k (t) = U k (t)dt + axdiW^t, dU k (t) = a 2 d{Wl) t) k = 1, n, 
X(0) = x, U(0)=u, t>0, 

where X(t) and U(t) are considered as random variables with given 
initial distributions, (X(t),U(t)) runs in the phase space R n x R™, o\ 
and 0"2 are nonnegative constants such that |a| ^ (o = (01,02)) 
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and (W l ) t = (W{, . . . , Wfyt, z = 1, 2, are independent n - dimensional 
Brownian motions. 

Let P(t, x, u) be a probability density of joint distribution of random 
values (X, U), subject to initial data 

n 

(1.3) P (x, u) = 5(u- u (x))po(x) = Y[ S(u k - (u {x)) k )p (x), 

k=l 

where po(x) is a bounded nonnegative function from L} oc {R n ). 
Let us introduce functions 

f uP(t,x,u)du 

(1.4) p a (t,x) — I P{t,x,u)du, u a (t,x) 



J P(t,x,u)du 

IK" Rn 

Notice that u a (0, x) = uo(x). Certain properties of u a (t,x) have been 
established in [2] and [3]. 

The density P = P(t, x, u) obeys the Fokker-Planck equation 



dP 



" d « i a 2 Ai 2 a 2 



p. 



subject to initial data (II. 3p . 

We apply the Fourier transform with respect to variables x and u in 
(11. 5p . (II. 3p and obtain the following Cauchy problem for the Fourier 
transform P = P(t, A, £) of P(t, x, u)\ 

(1-6) ^ = -^ (T 'l A l 2 + ^l^ i3+ ( A '^' 

(1.7) P(0,A,O = J e- l(x ' s) e~ i ^ Ms)) p (s)ds, A^6l". 

R n 

Equation (11.61) can easily be integrated and we obtain the solution: 
(1.8) 

P(t, A, = P(0, A, e+At) exp {-\al\\\H - ^ (| \\H 2 + 3f (A, + 3|£| 2 )) 

The inverse Fourier transform (in the distributional sense) allows to 
find the density P(t,x,u): 

P(t,x,u) = j±s J je^e^P(t,X,Od\d^ = 



(2k) 



2 n. 
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f r \u (s)-u\ 2 6\(u {s) + u)^+s-x\ 2 

Po(s) / / e^V^e 12CT i t+CT 2 2 ' 3 d\d£ds 



\u Q (s)—u\ 2 6|(uq(s)+u)^ + s — x\ 2 



(1.9) =C / p (s)e 2t °i 12 ^ t+CT 2* 3 tZs, 



where 

A = A(t, A, s) = 1 A+^= 

2V3 \JV2aft + er^ ^ 



S = S(f, C, s) = (7 2 Vte+^^ f y At 2 + zmq(s) - , C 
From (11. 9p we have 



na 2 t a/12ct 2 + cxft 2 



AT \ f 3\u (s)t + s-x\^ 



(1.10) p CT (t,x) = — = / p (s)e 2 «H+- 2 2 < 2 ) rfs 

1 (v/27rf(3a 2 + a 2 2 t2) 1 ' 



(1.11) u CT (t,x) 



3|u{)(s)t + s — x| 

/ F(t,x,s)p (s)e~ 2 '(^+- 2 2 * 2 ) ds 

1 R™ 



/ p (s)e 2t < 3 ^ +CT 2 2t2 ) rf s 

where F(t, x, s) = 6afuo(s) — crft(ttto(s) + 3(s — x)). 

2. Properties of the "perturbed velocity" 

Proposition 2.1. Let uq and po > be functions of class C 1 (M n ) (1 
Cb(R n ). Ift*(uo) > is a moment of time such that the solution to the 
Cauchy problem with the initial condition uq keeps this smooth- 

ness for < t < £*(tto) < +oo, then u a (t,x) tends to a solution of 
problem U.l\) as \a\ — > for any fixed (t,x) G R n+1 , < t < i*(u ). 

Proof. Let us denote by J(mo(x)) the Jacobian matrix of the map x i — )• 
Mo(x). As it was shown in [I] (Theorem 1), if J(u (x)) has at least 
one eigenvalue which is negative for a certain point x 6 M n , then the 
classical solution to (II. ip fails to exist beyond a positive time t*(uo). 
Otherwise, i*(tto) — °°- The matrix C(t, x) = (/ + tJ(wo(x))), where 
/ is the identity matrix, fails to be invertible for t = U(uq). 
From (II .lip and (11.101) we have 

u a {t,x)p a (t,x) = (u*(*,a:) -u"(t,x))p a (t,x), 
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where 



ui(t,x) = K / u (s)p (s)e ™&°i+°W) ds 



<\t, x) = I (tuo(s) + s - x) Po ( S y da, 

2(3a( + a$t z ) J 

K = — , Using the weak convergence of measures and 

the fact that po and uq are continuous and bounded, we get for u l a it, x): 

(\ n 3\uQ(s)t+s — x\^ 

lim = ^ — al „^.^._ T|2 - 

H->o 



(\ « 3|ii (s)t + s-a:| 2 



/ u (s)p (s)5 p (t tX>s )ds 

/ Po{s)8 P (t,x,s)ds 

t" 

with p(£, x, s) = «o(s)t+s — x , where S y is the Dirac measure as y G M n . 
We can then on the basis on the invertibility of C (t, x) use locally the 
implicit function theorem and find s = St, x {p)- 
Therefore, 

/ Uo(s t , x (p))po(s tjX (p)) det(C(t, st^ip)))- 1 5 P (ds t!X ) 

liixi (tj = ~ = 

M-»° a ' J Po(s t ,x(p)) det(C(t, s t>x (p)))~ l 5 P (ds tjX ) 

= uo(s t>x (0)). 

Now we prove that lim u IJ (t, x) = 0. We change variable 5 = ^"°^ t+ = 

|cr|-s>0 ^/t(3(T 2 +<x 

and obtain 

/ cr|tp (s) exp (—s 2 /2)ds 
lim u^it.x) = lim ; : = 0. 

Let us introduce the new notation so(t, x) = s ttX (0). Then the fol- 
lowing vectorial equation holds: 

(2.1) u (s (t, x))t + s (t, x) — x — 0. 

The function u(t, x) = u (s (t,x)) satisfies the Burgers equation. To 
check this fact it is enough to differentiate (12. ip with respect to t and 
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Xj and find so,t an d s o,x- The more detailed proof you can find in pQ in 
Proposition 1. □ 

It is important to note that So(t,x) is unique for all t for which the 
solution to the Burgers equation u(t, x) is smooth. 

Remark 2.2. Proposition 12.11 can naturally be extended to the class of 
functions p such that there exists a sequence p e e C 1 (lR n ) fl C^iW 1 ) 
converging to p as E — > almost everywhere. In this case u a (t,x) 
tends to a solution of problem (II. ip as a — » almost everywhere on 
G M n+1 , < t < U{uq). 

Theorem 2.3. The scalar function p a (t,x) and the vector-function 
u a (t,x), defined in \1.4\ ), solve the following system: 

(2.2) — + div x {p a u a ) = -a l ^ — T , 

k=i k 



d(p a u a>i ) 1 2 f 

(2.3) — h div x {p a u aii u a ) = ^ a i2^ 



d 2 {p a u (7 ,i)) 



, .....,..,.„...,,..,„ go 

k=l ft 

where i — 1, .., n, t > and I a = J (ui — u at i)((u — u a ), V x P(t, x, u))du. 



Proof. Integrating ( 11. 5p with respect of u we get: 

dp a 1 2 A3 2 p ff 1 2 f fd 2 P(t,x,u) 

jU_1 fe t. — 1 ft 

The last integral is equal to zero, because P(t,x,u) and all its deriva- 
tives tend to zero as \u\ — > oo. Thus, we obtain (12. 2p . 

To prove (12.31) we note that the definitions of u a (t,x) and p a (t,x) 
imply 



d(p a u a ,i) d 



dt 



— J UiP du = J UiP t du = —J Ui(u,V x P)du+ 



' 1 d 2 P . 1 „ A /• d 2 P , 

w 7 - - „ aw, 



(2.4) + -a? £ y d« + -a 2 2 £ y 



ft — J-TED-n. ft — J- TO) 77. 



where P$ = 4rP. Early we noted that the last integral is equal to zero. 



dt J 

Further, for % = 1, ... n we have 



div x (pu aii u a ) = u ati div x (p a u a ) + u a d x (p a u aii ) - u (T u (T)i d x p 
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(2.5) = J u a . ii (u,V x P)du + J Ui{u a) V x P)du- J u aji (u a ,V x P) du, 

R n R n R n 

where i, k = 1, n. Equation f l 2 . 3 H follows immediately from ( 12.4ft and 
( 123|) . □ 

Let us set p(t,x) = lim p a (t,x) and u(t,x) = lim u a (t,x). 

\<r\-*0 |cr|— >-0 

Theorem 2.4. Assume that (p(t,x),u(t,x)), the limits of (p a ,u a ) as 
\o~\ —7- 0, areC 1 -smooth bounded functions for (t, x) G := [0, £*(«())) X 
IR n , t*(«o) < oo. Then they solve in Q the pressureless gas dynamics 
system 

(2.6) d t p + div x (pu) = 0, d t (pu) + V x (pu ® u) — 0. 

Proof. As follows from Proposition 12.11 the function -u(i, x) is a C 1 - 
solution of the non-viscous Burgers equation. Further, (12. 2 p is a linear 
parabolic equation with respect to p a , hence the limit as | o \ — >■ re- 
duces it to the first equation in ( 12 .6p (continuity equation). The second 
equation in (I2.6P is a corollary of the non- viscous Burgers equation and 
continuity equation for smooth solutions. □ 

Remark 2.5. Proposition 12 .41 implies that the integral term on the right- 
hand side of (12. 3p vanishes as \o\ — > in the case of smooth limit 
functions p and u. 

3. A REGULARIZATION OF THE PRESSURELESS GAS DYNAMICS 

After the moment of shock formation in the solution to the Burgers 
equation the limit of (I2.2p - (l2.3p as |a| — > is not equivalent to the 
pressureless gas dynamics system ( 12. 6p . the integral term in (12. 3p does 
not vanish and can be interpreted as a gradient of special pressure (spu- 
rious pressure). It is well known that for the pressureless gas dynamic 
system after the moment of singularity formation there arise a question 
how to define a generalized solution, and the answer depends on the 
prescribed type of interaction between particles. For the sticky par- 
ticles model, where two particles move together after interaction, one 
need to introduce a strongly singular solution with a delta-singularity 
in the density component. If we assume that particles do not interact, 
we get the generalized solution in the sense of free particles, whose sin- 
gularities are only shocks. In this sense we can call the limit of system 
( I2.2p - (l2.3p as \a\ — > a regularization of the pressureless gas dynamics 
equations. There exist an algorithm how to pass from the model of free 
particles to the model of sticky particles (see [1] for details in the case 

0- 2 = 0). 
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For discontinuous initial data analogously to [T] we introduce a gen- 
eralized solution in the sense of free particles: 

Definition 3.1. We call the couple of functions (ppp(t, x),Upp(t, x)) 
the generalized solution to the Cauchy problem for (I2.6P in the sense 
of free particles (FP-generalized solution) subject to initial data 
(po(x),u (x)) e Lf oc (M n ) nL 00 (R n ), if for almost all (t,x) GK + xl B 

PFp(t,x) = limflim p £ a (t, x)), UFp(t,x) = limflim u £ (t, x)), 

£— >0 a— >0 £— >0 a— 5>0 

where (p^it, x), u £ a (t, x)) correspond to initial data p e Q = r] £ *p , u £ Q = 
7] £ *uq , where r) £ (x) is the standard averaging kernel. 

In [1] for o 2 = we solve the Riemann problem in the sense of free 
particles, this can be done analogously in the general case a 2 > 0. 

It is very important that the double limit procedure used in the 
definition of the FP-solution plays role only for the case of a central 
rarefaction wave, where it helps to find a unique stable solution in the 
vacuum domain (we refer to [TJ again). For the compression wave the 
FP-solution can be computed directly by formulae (11.101) . (11.111) using 
only one limit pass \a\ — > 0. Thus, after the moment of the shock 
formation t*(«o) in the Burgers equation the solution to the respective 
conservation law 
(3.1) 

d t p + div x (pu) = 0, d t (pu) + V ' x (pu <g> u) + Vp = 0, Vp = lim l a 

H->o 

can be found from (11.101) . (II. lip as |cr| — > 0. By definition the solution 
to (13.11) is the FP-solution to (12.61) . thus (13.11) can be considered as a 
regularization of (12.61) . 

4. The ID model for polytropic gas 

Let us assume that the gas that obeys system (13.11) is polytropic, 
i.e. the state equation is p = Ap 1 , 7 is the adiabatic exponent, A is a 
positive constant. Therefore we have to supplement (13. ip with equation 

d t p + ud x p + •ypdxU = 0. 

We consider the ID Riemann problem in the compression case U2 < 0: 
po(^) = pi + p2@(x), u (x) — U\ + u 2 0(x), where p i; U{, i = 1,2 are 
constants, In pQ the solution was obtained (for 02 = 0, the general case 
is analogous): 

p(t, x) = p x + (pi + p 2 )9(x - uit) - pi9(x - (ui + u 2 )t), 

u(t, x) — ui + f 1 + P2 u 2 6(x - u\t) - - Pl — u 2 9(x - (ui + u 2 )t), 
2pi + p 2 2pi + p 2 



s 
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p (t, x ) = pl } pl + p2 K 2 2 (9(x - { Ul + u 2 )t) - 6{x - u x t)). 

2pi + p2 

Proposition 4.1. In assumption that the state equation is polytropic, 

the adiabatic exponent 7 — > 3 as the ratio — — > 0. 

P2 

Proof. We have three Hugoniot conditions, that follow from the mass 
conservation, momentum conservation ( f)3.ip ) and the conservation of 

1 p 

total energy E = -pu 2 H . We denote by /_ and /+ the values of 

2 7 — 1 

f(x) before and after shock, respectively. For the shock moving with 

the speed D = u\ + u<i we have p + = 2p\ + p 2 , p- = pi, u + = u\ + 

Pi + P2 (h{fh +P2) 2 n ™ ... ,„ • , 

U2, w_ = Ui, p + = ui, P- = U. the third Hugoniot 

2pi + p 2 2pi + p 2 

condition gives D = - — M an( j ^ — ^El ^ 2 _ This implies the 

[S] Pi + P2 

proposition. The consideration of the second shock gives the same 
result. □ 

The adiabatic exponent 7 = 1 + -, where n is the space dimension, 
corresponds to the monoatomic gas. 
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